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ABSTRACT 

The type Ic supernova (SN) 2002ap is an interesting event with very broad spectral features like the famous 
energetic SN 1998bw associated with a gamma-ray burst (GRB) 980425. Here we examine the jet hypothesis 
from SN 2002ap recently proposed based on the redshifted polarized continuum found in a spectropolarimetric 
observation. We show that jets should be moving at about 0.23c to a direction roughly perpendicular to us, and the 
degree of polarization requires a jet kinetic energy of at least 5 x 10 50 erg, a similar energy scale to the GRB jets. 
The weak radio emission from SN 2002ap has been used to argue against the jet hypothesis, but we argue that 
this is not a problem because the jet is expected to be freely expanding and unshocked. However, the jet cannot 
be kept ionized because of adiabatic cooling without external photoionization or heating source. We explored 
various ionization possibilities, and found that only the radioactivity of 56 Ni is a plausible source, indicating that 
the jet is formed and ejected from central region of the core collapse, not from outer envelope of the exploding 
star. Then we point out that, if the jet hypothesis is true, the jet will eventually sweep up enough interstellar 
medium and generate shocks in a few to 10 years, producing strong radio emission that can be spatially resolved, 
giving us a clear test for the jet hypothesis. Discussions are also given on what the jet would imply for the GRB-SN 
connection, when it is confirmed. We suggest existence of two distinct classes of GRBs from similar core-collapse 
events but by completely different mechanisms. Cosmologically distant GRBs having energy scale of ~ 10 50-51 
erg are collimated jets generated by the central activity of core collapses, associated with 56 Ni ejection along 
with the jets. SN 2002ap can be considered as a failed GRB of this type with large baryon contamination. On 
the other hand, much less energetic ones including GRB 980425 are rather isotropic, which may be produced by 
hydrodynamical shock acceleration at the outer envelope. We propose that the radioactive ionization for the SN 
2002ap jet may give a new explanation also for the X-ray line features often observed in GRB afterglows. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — polarization — stars: circumstellar matter — 
supernovae: individual (SN 1998bw, SN 2002ap) 



1. INTRODUCTION 

Type Ic supernova (SN) 2002ap has attracted particular atten- 
tion since its discovery by Y. Hirose in January 2002, because 
of its relatively close distance (about D = 7.3 Mpc, Sharina, 
Karachentsev, & Tikhonov 1996; Sohn & Davidge 1996) and 
its broad-line spectral features (Kinugasa et al. 2002) that are 
considered as a signature of very energetic supernovae. Such a 
supernova population, often called as hypernovae whose proto- 
type is the famous type Ic SN 1998bw, have explosion energy 
more than 10 times larger than the standard energy (~ 10 51 erg) 
when spherical symmetry is assumed (Iwamoto et al. 1998; 
Woosley et al. 1999; see also Hoflich, Wheeler, & Wang 1999 
for asymmetric modeling of these events with less extreme ex- 
plosion energies). The apparent association of a gamma-ray 
burst (GRB) 980425 with SN 1998bw makes these mysteri- 
ous events even more interesting in the context of the possible 
SN/GRB connection. 3 Mazzali et al. (2002) presented photo- 
metric and spectroscopic modeling of SN 2002ap assuming a 
spherical explosion, and indicated that the explosion occurred 
at Jan 28±0.5 UT, with a kinetic energy of about (4-10) x 10 51 
erg and the progenitor is a C+O star whose main sequence mass 
is ~ 2O-25M . It seems that an interacting binary is more 
likely for a star of this mass scale to lose its hydrogen and he- 



lium envelope, but theoretical and metallicity uncertainties do 
not reject a single Wolf-Rayet (WR) star as another possible 
progenitor (Smartt et al. 2002). 

In contrast to SN 1998bw / GRB980425, SN 2002ap was not 
associated with a GRB to the sensitivity of IPN (Hurley et al. 
2002; but see also Gal- Yam et al. 2002). On the other hand, 
spectropolarimetric observations of SN 2002ap (Kawabata et 
al. 2002; Leonard et al. 2002; Wang et al. 2002) give an inter- 
esting hint for hidden energetic ejecta. Kawabata et al. (2002) 
noticed that the spectral shape of polarized continuum observed 
by Subaru around 10 Feb (i.e., ^13 days after the explosion) ap- 
parently looks like the original unpolarized spectrum, but red- 
shifted by z = 0.3 (A re dshifted/A = 1 +z) and the ratio of the polar- 
ized to unpolarized flux is fp = 0.0018 (in f\). The polarization 
angle (PA) is different from line polarization at this epoch or PA 
in their second observation in March (40 days after the explo- 
sion). Interestingly, they got a consistent PA and wavelength- 
independent polarization from February to March, that can be 
explained simply by an asymmetric photosphere as often seen 
in supernova spectra, after they subtracted the redshifted polar- 
ized continuum in February observation. If it is not a chance 
coincidence, this result can be explained by an asymmetric su- 
pernova photosphere and a jet moving at a much higher speed 
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3 GRB 980425 was exceptionally faint one compared with those found at cosmological distances. However, after the submission of this paper, a supernova SN 2003dh 
having similar spectrum to SN 1998bw was discovered to be associated with a nearby GRB 030329 having normal luminosity (Stanek et al. 2003; Hjorth et al. 2003), 
finally confirming the connection between GRBs and energetic supernovae. 
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(~ cz ~ 0.3c) than the supernova photosphere (Kawabata et al. 
2002). Following this suggestion, Leonard et al. (2002) con- 
firmed the resemblance between polarized and redshifted spec- 
trum, by an independent data taken by Keck, although statistical 
significance of this resemblance is difficult to assess. 

As we will show in §2, jet mass of Mj et ~ 0.01M Q and jet 
kinetic energy of £j et <~ 10 50-51 erg are required to explain the 
degree of polarization. It is quite interesting in the perspective 
of GRB-SN connection, to note that the inferred jet energy is 
very close to the standard energy scale of collimated jets sug- 
gested for GRBs (Frail et al. 2001; Panaitescu & Kumar 2001). 
It is well known that an ultrarelativistic outflow with a Lorentz 
factor 7 > 100-1000 is required for successful GRBs to avoid 
the compactness problem (Goodman 1986; Paczyhski 1986). It 
may be achieved by production of a fireball, where enormous 
energy is injected into a clean region with very few baryon con- 
tamination [~6x lO- 6 (£/lO 51 erg)( 7 /lOO)- 1 M ]. Then it is 
naturally expected that there may be events with similar jet en- 
ergy but with much larger baryon contamination, and hence low 
expansion velocity and no gamma-ray emission, often called 
as "failed GRBs" or "dirty fireballs" (e.g., Huang, Dai, & Lu 
2002). The energetic jet inferred for SN 2002ap may be the 
first detection of a failed GRB. On the other hand, even suc- 
cessful GRBs may also have jets with low velocity and high 
baryon load, which are produced in the process of jet acceler- 
ation in addition to the ultra-relativistic component responsible 
for GRBs. SN 2002ap may have been a successful GRB, but 
the jet was not directed to us. 

However, this jet hypothesis has been questioned at a few 
points. Radio emission is thought as an indicator for the exis- 
tence of fast moving ejecta, since it would produce nonthermal 
synchrotron emission by interaction with circumstellar mat- 
ter (CSM) or interstellar matter (ISM). However, radio emis- 
sion of SN 2002ap is much (by more than 3 orders of mag- 
nitude) weaker than that of SN 1998bw, and Berger, Kulka- 
rni, & Chevalier (2002, hereafter BKC02) have shown that 
the observed radio emission can be explained by synchrotron 
radiation produced by a spherical ejecta expanding at ~ 0.3c 
with a total energy of nonthermal electrons of E e <~ 1.5 x 10 45 
erg, assuming energy equipartition between electrons and mag- 
netic fields. These numbers should be compared with those 
for SN 1998bw, i.e., relativistic shock speed with a Lorentz 
factor 7 - a few and E e ~ 10 49 erg (Kulkarni et al. 1998). 
Based on this result, BKC02 argued that the energetic jet with 
£jet ~ 10 51 erg proposed by Kawabata et al. (2002) should have 
produced much stronger radio emission than observed. Wang 
et al. (2002) presented VLT spectropolarimetric observations 
including earlier epochs than Keck and Subaru, and found that 
the continuum polarization evolved from nearly zero on 3 Feb 
to 0.2% on 10 Feb, which is contrary to what is expected from 
a simple jet model, since scattering efficiency should be higher 
in earlier epochs when jet location is closer to the star. Finally, 
it is uncertain whether the jet material is kept highly ionized in 
spite of the expected rapid adiabatic cooling. 

The main purpose of this work is to examine whether the jet 
hypothesis is physically tenable especially against the possible 
difficulties mentioned above. We found that it is in fact physi- 
cally possible that the jet exists but has been spirited away from 
intensive observational efforts made so far, except for the Sub- 
aru spectropolarimetry. However, it is still a hypothesis, and 
we need a further observational test to prove this. Fortunately, 
we show that there is a good test for this hypothesis by future 



observation; we point out that a long-term radio monitoring of 
this object should find re-emergence of radio emission in a few 
to 10 years, for which the jet expansion should be easily re- 
solved by VLBI imaging. Then, we will give discussions on 
what would be the implications for the GRB-SN connection if 
the jet is confirmed in the future. Here we also propose a new 
possible explanation for X-ray line features often observed in 
GRB afterglows, which is inspired by the results of the paper. 

The paper is organized as follows. We give an estimate of 
the jet mass and energy with a fully relativistic treatment in §2. 
We discuss the physical condition of the jet including ionization 
sources in §3, and give a detailed modeling of radio emission 
from the interaction between the jet and CSM/ISM in §4. §5 
is for discussion on the GRB-SN connection and X-ray lines in 
GRB afterglows, and summary and conclusions are presented 
in §6. 

2. JET MASS AND ENERGY ESTIMATION 

We assume that the jet is sufficiently collimated and hence 
we can define a single scattering angle of photons scattered by 
the jet, o bs, which is the same with the direction angle of the 
observer measured from the jet direction. We also assume that 
the jet is optically thin for the electron scattering, which will be 
checked later. If the jet is optically thin, the jet mass and energy 
estimates in this section do not depend on the jet opening angle. 
We use a notation that x and x' are the quantities in the rest- 
frames of the supernova/observer and the jet, respectively. [The 
heliocentric redshift of the host galaxy is +631km/s (Smartt & 
Meikle 2002) and can be neglected.] According to the standard 
Lorentz transformation, the energy of original photons from the 
supernova has a relation 

e org = 7(l-/3)eorg , (1) 

and for the energy of scattered photons by the jet: 

e sc = 7( 1 - /3 cos 6> obs ) e sc , (2) 

where (3 and 7 are the bulk velocity and Lorentz factor of the 
jet, respectively. From these relations, and using e' OIg = e' sc for 
the Thomson scattering, we have the relation between the orig- 
inal and scattered photon energies: 

esc = 1 = 1-/3 



(3) 



corg 1+Z l-/3cOs6> obs ' 

The inverse Lorentz transformation of the equation (2) leads 
to: 



e S c = 7(l + /3cos6£ bs ) 4 , 
and hence 6* bs and O' obs are related as: 



(4) 



7 2 (l-/?cos0 obs )(l + /3cos0 obs )=l • 



(5) 

We should consider three time scales: the time when we ob- 
serve a scattered photon (f b s , measured from the arrival time 
of unscattered photons emitted at the core-collapse date), that 
when the photon is scattered by the jet (f sc ), and that when the 
photon is originally radiated from the supernova photosphere 
(fad). The last two are measured from the core-collapse date in 
the supernova/observer restframe. From a geometrical calcula- 
tion, we find 



frad — 



1 



l-/3cos# obs 



fobs 



1 



1-/3 

l-/}COS6> bs 1+Z 



(6) 
(7) 
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Let Forg&c) = £or g (W)/(47rr? t ) be the original flux from the 
supernova at the jet location, rj et = c(3t %c , and the flux in the jet 
restframe is 

F ' rg = 7 2 (l-/3) 2 J F org . (8) 

The luminosity of scattered light per unit solid angle is given 
as: 

_ N do(6' ohs ) 

dn - e dn org ' ( > 

where N e = Mj et / e i / (p e m P ) is the free electron number in the jet, 
/Lt c the nucleon to electron number ratio, m p the proton mass, 
and /ei the fraction of free ionized electrons. Assuming that the 
jet material is mostly heavy element, e.g., C+O, we set \i e = 2. 
Here, da(9' oh ^/dVt = (3/167r)ar(l + cos 2 # obs ) is the cross sec- 
tion of Thomson scattering for unpolarized light. The scattered 
luminosity dL' sc /dft is related to that in the supernova/observer 
restframe as (e.g., Rybicki & Lightman 1979): 

dn 7 4 ( 1-/3 cos <9 obs ) 4 dn } 
The ratio of polarized to unpolarized flux, f P , is given as: 

where II = (1 - cos 2 # obs )/(l +cos 2 8' obs ) is the degree of polariza- 
tion of scattered wave, and the last factor of (1 + z) _1 is coming 
from the definition of f P by the ratio of flux per unit wavelength, 
fx- 

From eqs. (6) and (7), the fractional time delay of scattered 
photons from direct unscattered photons is (f bs -W)Aobs = 
1 — (1 +z) _1 = 0.23, which is small and independent of unknown 
#obs- Since the luminosity and spectrum of the supernova are 
not expected to change significantly within these time scales, 
we expect that the scattered spectrum is similar to the unscat- 
tered one except for the redshift, as observed. Therefore we do 
not have to take into account the luminosity and spectral evolu- 
tion of supernova, i.e., t m & ~ f Q b s . In the top and middle panels 
of Fig. 1, we show (5, 7, Mj et , and £j et = Mj et c 2 (j- 1) required 
to reproduce the observed redshift (z = 0.3) and degree of polar- 
ization (f P = 1.8 x 10" 3 ), as a function of the jet viewing angle, 
#ob s - Here, we have assumed f e \ = 0.3 and f b s = lOd, and scaling 
of the results by different values of these parameters is obvious. 

It is likely that there is another jet in the opposite direction 
from the supernova, in addition to the jet considered so far, and 
the opposite jet should also produce another redshifted polar- 
ization component. Since the data show only one redshifted 
component, the contribution from this opposite jet must be with 
similar polarization degree and redshift to the original jet, or 
negligibly small due to too small fp or very large redshift. In the 
bottom panel of Fig. 1, we plot the redshift (z op ) and the polar- 
ization degree (/p )0 p) by the opposite jet having the same jet ve- 
locity and mass but # b s , p = 180°-6'obs- This gives a constraint 
°f #obs ^ 100°, otherwise we should have observed another po- 
larized continuum component with larger fp and smaller red- 
shift than the observed ones. 

The jet velocity is roughly constant at (3 ~ 0.2 for 9 \, s > 90°, 
but it becomes more relativistic with decreasing 9 bs at o b s < 
90°, because the redshift effect of the jet motion is compen- 
sated by the blueshift to the observer. The jet mass and energy 
rapidly becomes larger with decreasing 9 D \, S , mainly because 
of less efficient polarization and larger rj et oc f sc 3> fobs- Small 
#obs(;5 60°) seems not favored from energetics, since it requires 



jet energy of more than 10 erg. These considerations lead to 
a conclusion that the jet directions must be close to 9 \, s ~ 90° 
(probably within ± 10-20°) and both the two jets contributed 
roughly equally to the observed polarization. Therefore we as- 
sume two jets with # b s = 90° in this work. Then we found 
(3 = 1 - (1 + z) _1 = 0.23, and the jet mass is 

Mjet = 0.01l(^^V^)Y^0 M , (12) 
J V0-0018y VIOd/ \0.3J ' v ' 

where Mj et is re-defined as the mass of each jet. Therefore, ob- 
served redshifted polarization can be explained if the jet ma- 
terial is modestly ionized, with the kinetic jet energy £j et ~ 
5 x lO 5O (M jet /O.OlM )(/3/O.23) 2 erg. In the following of this 
paper we consistently use o b s = 90°, f sc = f b s , /3 = 0.23, and 
M jet = 0.01M Q . 

3. PHYSICAL CONDITIONS OF THE JET AND IONIZATION 

3.1. Jet is freely expanding 

First we consider the fact that the radio emission from SN 
2002ap was very weak. If a considerable part of the jet kinetic 
energy was converted into nonthermal electrons via shock ac- 
celeration, inevitably there must be very strong radio emission 
which should have been even stronger than SN 1998bw. How- 
ever, if the amount of CSM swept-up the jet is much smaller 
than the jet mass, the jet feels almost no deceleration and the 
majority of the jet material remains unshocked. We do not ex- 
pect radio emission from such an almost freely expanding jet, 
and we only expect radio emission by CSM swept-up by the 
jet. This emission and total energy of radio emitting electrons 
(<~ 10 45 erg) are simply related by the jet velocity and CSM 
density, and not related with the total jet mass and kinetic en- 
ergy. Stellar wind mass-loss rate of Wolf-Rayet stars, which 
are considered as a possible candidate of the SN Ic progeni- 
tors, is typically M w ~ 10~ 6 -10~ 5 MQ/yr with a wind velocity 
of V w ~ 10 3 km/s (McCray 1983; Garcfa-Segura, Mac Low, & 
Langer 1996; Garcfa-Segura, Langer, & Mac Low 1996). It 
is generally assumed that the CSM around radio supernovae 
has a stellar wind profile, i.e., pcswiir) = M w /(4-irr 2 V w ). Then 
the swept-up mass by the jet becomes M sw = bM w r ]et /V w = 
1.9 x lO- 7 fc-iM w -6V-^i O M , where b = 0.1Z?_i = fi je t/47r is 
the beaming factor of the jet opening angle, t\o = f o b s /10d, 
M w _ 6 = M w /(lO- 6 M /yr), and V W)3 = V lv /(10 3 km/s). This is 
much smaller than the jet mass and hence the jet is not deceler- 
ated. 

Radio emission from SN 1998bw and 2002ap can be ex- 
plained by isotropic high speed ejecta interacting with the CSM 
density consistent with the wind parameters similar to the above 
values (Kulkarni et al. 1998; Li & Chevalier 1999; BKC02). 
[The speed of ejecta is, however, considerably different for 
these two; 7 e j ~ a few for the former but u e j ~ 0.3c for the lat- 
ter.] Since the inferred jet velocity of SN 2002ap is close to v e j, 
and the swept-up CSM mass by the jet should be smaller than 
that by isotropic ejecta because of the collimation, it looks rea- 
sonable that the radio emission from the external shock front of 
the jet is equal to or smaller than the observed radio emission. 
We will present more detailed radio emission modeling in §4. 

Therefore the weak radio flux from SN 2002ap does not im- 
mediately exclude the jet hypothesis, if it is expanding freely. 
It may also be useful to recall that the kinetic energy of su- 
pernovae (~ 10 51 erg) is hardly converted into radiation, but 
supernovae are heated and shining by radioactivity. However, 
free expansion raises another problem because of the expected 
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rapid adiabatic cooling. The jet material must be ionized at 
least modestly, but as we will show below, the adiabatic cool- 
ing is so strong that the temperature of the jet material is likely 
lower than that necessary to keep it ionized. Therefore we need 
an external radiation or heating source of ionization. We will 
discuss this issue in detail below. 

3.2. Optical depth and thinning burst 

First we check the optical depth of the jet material to the 
photon-electron scattering. The jet must be mostly transparent 
to radiation at the day ~ 13, for the jet material to contribute 
the scattering and polarization efficiently. It can be written as: 



/elMjet <?T nnC o f . -1 ,-2 

7jet ~ . , , = 0.088 /el b_ x t m 

fj, e m p ^br L 



(13) 



Therefore the jet becomes optically thin at a time f t h ~ 

1/2 -1/2 

3.0 / el o_j days. This thinning time is close to the epoch 
when the continuum polarization evolved from nearly zero to 
the 0.2% level in the VLT observation, and hence indicating 
that the initially unobserved polarization can be explained by 
high optical depth of the jet. 

In analogy to the fireball theory for GRBs, we expect a burst 
of radiation when the jet becomes optically thin (Meszaros, La- 
guna, & Rees 1993), and here we examine how much radiation 
we expect from this. Suppose that the jet is in thermal equi- 
llibrium and the total internal energy is comparable with the 
kinetic energy of the jet, £j et , at the initial radius r,. If the jet 
interacts with stellar envelope and dissipate its kinetic energy, 
the radius of the progenitor C+O star, r, ~ 10 10 cm is a rea- 
sonable choice, while another possibility is r, <~ 10 6 cm if the 
jet is directly emitted from the central compact neutron star or 
black hole. The initial temperature 7} is determined by the in- 
ternal energy density t/, = Ej et /Vi, and both the internal energy 
and temperature adiabatically decrease as oc (rjet/V,) -1 , where 
Vi = A-xQr] is the initial volume of the jet material and ( is the 
fractional thickness of the jet shell. Here we assumed that the 
jet is homologously expanding like supernova ejecta, and hence 
V oc r? t . Then we find the internal energy and temperature at 
the thinning time fa as: 

E th ~ 2.8 x 10 45 /ei V yf r M0 erg 



(14) 
(15) 

where C-i = C/0-1 and r L1Q = r,/(10 10 cm). The flux 
of this thinning radiation is ~ E t h/(4-nD 2 fa) = 1.7 x 
10~ 12 /d 1 ^-i r i,io erg cm" 2 s _1 , which is sufficiently smaller than 
the observed UV-optical-IR flux of ~ 2 x 10" 10 erg cm" 2 s _1 at 
the day 5 (Sutaria et al. 2003). 

If there is no external heating or ionizing radiation source, 
the ionization balance of the jet would be determined by col- 
lisional ionization coefficient q co i by thermal electrons and ra- 
diative recombination coefficient a KC . Comparing coefficients 
of these processes (Lotz 1967; Nahar & Pradhan 1997; Nahar 
1999), the temperature must be higher than ~ 5 eV for the C+O 
matter to be ionized doubly or more. Considering the tempera- 
ture derived above, and the weak dependence on the unknown 
parameters, it seems unlikely that the jet is sufficiently hot to 
keep itself ionized. After the jet becomes optically thin, radia- 
tive cooling may further decrease the temperature. Then we 
need an external source of heating or ionizing photons. 

3.3. Photoionization 



To begin with, we estimate the total ionizing photon lumi- 
nosity required to keep all the jet material ionized (a similar 
argument used to derive the Stromgren radius). This is given 
by the recombination rate as: 

iph.rec ~ C^recWeMon 

= 5.6 x lO^llGVliTo 

/ « rec \-i fu_ 

llO" 11 cmV 1 / V 0.3 



(16) 



where n e is the free electron number density, N lon the number 
of ions, and yU; on = is the mean molecular weight of ions. 

Correspondingly, the ionizing flux in the jet must be stronger 
than fi on ~ L p h i rec/(47rfcr 2 et ), and the abundance ratio between 



species X ' + e 
N{X +i+l ) 



<-> X' is given as: 

f^ion^ion 



(^rectlp. 
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5 x 10" 18 cm 2 



(17) 



(18) 



Here, we used typical values of a mc and ionization cross sec- 
tion a loa for C + C +2 or + <-> +2 at temperature ~ 1 eV 
(Osterbrock 1989). Therefore, a necessary and sufficient con- 
dition for ionization is that the jet matter is radiated with a total 
photon luminosity given in eq. (16). 

First we consider the possibility of ionization by radiation 
from the supernova. The effective temperature inferred from 
optical colors of SN 2002ap at day <~ 10 is about T e ff ~ 6,000 
K (Gal-Yam et al. 2002), and the number of photons in the 
black-body tail above vr~6x 10 15 GHz, which is a threshold 
frequency of ionizing photons for CII and Oil, is quite small. 

There is an XMM-Newton observation of SN 2002ap on 
Feb 2 UT, i.e., about 5 days after the explosion (Soria & 
Kong 2002; Sutaria et al. 2003). The flux in 0.3-10 keV is 
1.07 x 10~ 14 erg cm" 2 s _1 , and the spectrum can be fit with a 
power-law with a photon index of a x = 2.6+qj (f v oc v~ ax+l ). 
The flux extrapolated down to UV bands may ionize the jet, and 
softer spectrum gives stronger ionization flux. 

However, following arguments give further constraints on the 
possible range of ax- Sutaria et al. (2003) ascribed this X- 
ray flux to be inverse-Compton scattering of the optical pho- 
tons, and in this case the extrapolated nonthermal flux down 
to the optical bands should be about TcsM^opt> where tcsm is 
optical depth of hot electrons in shocked CSM around the su- 
pernova and F opt the optical flux of the supernova (Fransson 
1982; Chevalier & Fransson 2001; Sutaria et al. 2003). The 
hot CSM gas responsible for X-ray emission is swept up ei- 
ther by the jet or by another isotropic fast ejecta as consid- 
ered by BKC02 (see also §4 for examination of these two 
possibilities from radio data). In either case, the velocity of 
the shock front is 0.23-0. 3c, and hence the location of X-ray 
emitting region is r x ~ 6.0 x 10 15 fio cm. Therefore we find 
tcsm = 2.0 x \Q~ 5 M w - b V~\t^bx, where bx is the sky cover- 
age of X-ray emitting region viewed from the supernova pho- 
tosphere; bx = 1 for the isotropic ejecta while bx = b for the 
jet. As we will show in §4, M w < lO~ 4 M /yr is required for 
successful modeling of the observed radio data. Then, compar- 
ing with the observed optical flux, ax < 2.6 and 2.2 is required 
for bx = 1 and bx = b = 0.1, respectively, even though softer in- 
dex is allowed by the observational error. BKC02, on the other 
hand, suggested that the X-ray flux is explained by the same 
synchrotron radiation as the radio observations. The radio and 
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X-ray fluxes at the day 5 are connected by a power law with 
a photon index of a' x = 1.5, and considering that the photon 
index changes by 0.5 below and above the cooling break fre- 
quency, the maximum photon index allowed in the X-ray band 
is ax < a' x + 0.5 = 2, irrespective of bx- 

Now we compare these constraints on ax with the range re- 
quired to ionize the jet. The observed X-ray flux at the day 5 
should not be much different at the day ^40-13, and we extrap- 
olate the X-ray luminosity down to the UV band and compare 
to Lph.rec- Note that only a fraction of the X-ray luminosity is 
directed to the jet material, and this fraction is given by ~ b/b x 
from a geometrical consideration. We found that the spectral 
index must be extremely soft as ax > 5 or 4 for bx = 1 or 
bx = b = 0.1, respectively, in order that the extrapolated flux 
down to vj is equal to Lph.rec- Therefore we can safely exclude 
the possibility that the nonthermal radiation producing the ob- 
served X-rays is ionizing the jet. 

Secondly we consider a possibility that a hot, UV-radiating 
star close to the SN 2002ap may ionize the jet. It is expected 
that SN 2002ap occurred in a massive star forming region where 
young massive stars are clustering. A close binary system is a 
candidate for the type Ic supernova progenitors (Nomoto, Fil- 
ippenko, & Shigeyama 1990), and it may provide even stronger 
ionization source. However, the total ionization luminosity 
given in eq. (16) is even larger by a factor of several than the 
ionization flux, ~ 10 49 5 s" 1 , above the frequency v T 10 15 
Hz for the most luminous and hottest stars (Schaere & de Koter 
1997). It should be noted that this luminosity is for all direction, 
but only the radiation within the solid angle of the jet viewed 
from the ionizing star is available for the jet ionization, which is 
expected to be a small fraction. If the region around SN 2002ap 
is filled up by radiation field with ~ F wn to a radius of rj et , the re- 
gion should have a luminosity of at least 47r,Fj on r? et , correspond- 
ing bolometric luminosity of 2.1 x lO 8 Z?Ij£l 1 1 /i74 f ro(/ei/O.3)L0, 
assuming the spectral energy distribution of the most luminous 
O stars. Such a huge luminosity is apparently ruled out by the 
prediscovery image of the SN 2002ap field reported by Smartt 
et al. (2002). To conclude, ionization by nearby young stars is 
impossible. 

3.4. Radioactive heating and collisional ionization 

Since photoionization of the jet seems difficult, the only way 
to ionize the jet is enhanced collisional ionization by external 
heating. If the jet is generated at the central compact object, it 
might include a significant amount of radioactive nuclei such 
as 56 Ni. Asymmetric explosion induced by the jet should also 
affect nucleosynthesis, and 56 Ni production along the jet direc- 
tion is enhanced (Nagataki 2000; Maeda et al. 2002). 56 Ni de- 
cays by electron capture and gamma-ray emission to 56 Co, with 
an exponential decay time scale of f N i = ?i /2 / ln(2) = 8.5 d and 
decay energy is eNi = 2.1 MeV. When the material is optically 
thick, the radioactive heat is quickly thermalized into optical 
radiation field, as generally seen for supernovae. On the other 
hand, if the material is mildly optically thin, the gamma-rays 
emitted by decaying 56 Ni scatter electrons with a probability 
~ Tjet, and since the gamma-ray energy is comparable with the 
electron rest mass, the scattered electrons acquire mildly rel- 
ativistic speed and energy. Such high energy electrons would 
lose their energy by ionization loss in the jet plasma, with a 
time scale of 

f;i ~ 2.0 x 10 4 (Vc) 4 4C-i£-i s, (19) 



where v e is the initial velocity of high energy electrons. Here 
we used the ionization loss formulae of Longair (1992) and the 
logarithmic factor is set to be 15. Therefore the energy de- 
posited by radioactive gamma-rays is used to ionize the jet ma- 
terial within the time scale of interest, giving an efficient ion- 
ization process. When optical depth is very low, this process 
would be dominated by positrons emitted from decay of 56 Co, 
which has a longer exponential lifetime of fc = 111.26 days 
and energy fraction given to positrons is 3.5% of the total de- 
cay energy (Arnett 1979; Woosley, Pinto, & Hartmann 1989). 

The ionizing balance is determined by the energy balance be- 
tween radioactive heating and recombination cooling (see also 
Graham 1988) as: 

/ Ni Mj e , exp(-f sc /f Ni ) 

rjet £Ni > a lec n e Ni on w , (20) 

^ Ni m p f N i 

where / Ni is the 56 Ni mass fraction in the jet, and w the ioniza- 
tion potential. The recombination rate coefficient depends on 
the electron gas temperature, which is determined by balance 
between the radioactive heating and cooling processes. We can 
estimate the minimum amount of 56 Ni by taking the minimum 
value of a rec as 

to > 0-68 tiCl*- 10 *'*^ ( 3xlo Xn3 g -0 (21) 

where the adopted value of a rec is minimum value of doubly 
ionized oxygen or carbon at temperature of ~ 10 4 K (Nahar 
& Pradhan 1997; Nahar 1999) and w 2 o = w/(20 eV). The jet 
may include significant amount of heavier nuclei that are dif- 
ficult to ionize for a fixed value of / e i, but on the other hand, 
it may also include considerable helium that is easier to ion- 
ize. The helium could be mixed from remaining helium layer 
of the progenitor, or it may be newly synthesized. Production 
of helium is also enhanced along the jet direction in energetic 
jet-like nucleosynthesis (Maeda et al. 2002). We also note that 
highly ionized heavy nuclei, such as 56 Ni, should produce ob- 
servable line emission in X-ray bands, and hence the observed 
weak X-ray flux gives a constraint on the species of ionized ele- 
ments. (See §5.2 for possible connection to X-ray line features 
often observed in GRB afterglows.) Whatever the jet composi- 
tion is, the above result indicates that, if the jet is kept ionized 
by radioactive heating, it must have a considerable amount of 
56 Ni (mass fraction of order unity). This estimate is, however, 
very uncertain especially about the composition of the jet, a rec 
and electron temperature. More sophisticated treatment is nec- 
essary to determine the ionization status, but it is beyond the 
scope of the paper. Therefore, it is difficult to conclude that the 
jet should be ionized, but it seems the best candidate of ioniza- 
tion process among others. 

The jet may be ionized by gamma-rays of 56 Ni decay leak- 
ing from the photosphere of SN 2002ap, even if the jet does 
not have radioactive nuclei. In fact, ionization of helium en- 
velope above the photosphere, which is required to explain the 
observed He lines in SN 1987 A and type lb supernovae, is as- 
cribed to the leaking gamma-rays from photosphere (Graham 
1988; Lucy 1991). The mass of 56 Ni produced by SN 2002ap 
is estimated to be 0.07 ± 0.02M Q from the light curve model- 
ing by Mazzali et al. (2002), which is larger than the jet mass. 
However, the efficiency for gamma-rays to hit the jet is reduced 
by the beaming factor, b, and it is further reduced by escaping 
fraction from photosphere. Although some supernovae, includ- 
ing SN 1998bw, showed evidence that a significant amount of 
gamma-rays are leaking in late phase (> 30 days) (Nakamura 
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et al. 2001; Patat et al. 2001), the leaking fraction should not be 
large in early phase of ~ 10 days, when the optical luminosity 
is still glowing up by diffusion of radioactive heat to the photo- 
sphere. According to the model of Mazzali et al. (2002), about 
10% of gamma-rays are leaking, most of which are produced by 
56 Ni outside the photosphere, at the day ~ 10 (K. Maeda & K. 
Nomoto 2003, a private communication). This should be con- 
sidered as an upper limit for the leaking fraction since the model 
assumes the maximally possible mixing, i.e., uniform distribu- 
tion of 56 Ni. [The best-fit model of Mazzali et al. (2002) has 
less significant mixing.] 

Therefore, leaking gamma-rays from or around the photo- 
sphere seem less efficient than 56 Ni in the jet, if the 56 Ni distri- 
bution is what is expected from isotropic modeling. However, 
if the explosion is very asymmetric due to the jet formation ac- 
tivity and considerable amount of 56 Ni is ejected outside the 
photosphere, it may ionize the jet. It should also be noted that 
SN 1998bw produced much larger amount of 56 Ni (~ 0.7M Q , 
Iwamoto et al. 1998) than SN 2002ap. We cannot reject that a 
comparable 56 Ni was produced also in SN 2002ap, but most of 
it is well outside the photosphere where optical depth is low and 
radioactive decay energy mostly escapes as gamma-rays, not in 
optical bands. Such 56 Ni could be missed in the modeling by 
Mazzali et al. (2002) based on optical observations. Such ex- 
treme mixing and distribution of 56 Ni is unlikely to occur sim- 
ply by hydrodynamical instability in C+O stars (K. Maeda & K. 
Nomoto 2003, a private communication, see also Shigeyama et 
al. 1990). Hence, significant ejection of 56 Ni from the stellar 
core by jet formation activity is again indicated. 

4. RADIO EMISSION BY SWEPT-UP MATERIAL 
4. 1 . Early emission by interaction with presupernova wind 

The observed radio emission is considered to be produced 
by shocked CSM swept-up by high velocity supernova ejecta. 
Then there are two possibilities: (i) the observed radio flux is 
generated by CSM swept up by the jet responsible for the red- 
shifted polarization, or (ii) the radio flux is from CSM swept 
up by isotropic supernova ejecta that is a different component 
from the jet, as considered by BKC02, and the radio emission 
by CSM swept up by the jet was weaker than observed. The 
former option predicts that the shock front of the radio emis- 
sion region is not decelerating, otherwise the shock generated 
in the jet material would overproduce much the observed ra- 
dio flux. On the other hand, since the expansion velocity of 
isotropic shell considered by BKC02 is similar to that of the 
jet, the mass of the isotropic ejecta in the latter option must be 
much smaller than the jet mass, otherwise the isotropic compo- 
nent would spoil the redshifted polarization produced by scat- 
tering in the jet. The mass of isotropic component can be as 
small as the swept-up CSM, that is much smaller than the jet 
mass as argued above, and it should be decelerated by swept-up 
CSM. BKC02 found that the radio data is not sufficient to con- 
strain whether the shock radius is decelerating or not, and hence 
cannot constrain this possibility. Here we examine the possibil- 
ity (i) by a detailed modeling of the observed radio emission 
with a collimated jet. 

Given the jet opening angle and the density of CSM (deter- 
mined by M w and V w ), we can calculate the mass and energy 
density of shocked CSM swept up by the jet moving at a con- 
stant velocity, 0.23c. We assumed the strong shock limit with a 
compression factor of 4, and temperature of the shocked CSM 
is calculated by the standard shock theory for supersonic pis- 



ton. Then we can calculate the synchrotron flux according to 
the standard formulae, if fractional energy densities of nonther- 
mal electrons (e e ) and magnetic field (e B ), electron power index 
(p,dN e /d-f e oc j~ p ), and the minimum Lorentz factor of non- 
thermal electrons (-y m ) are specified. We calculated the syn- 
chrotron flux taking into account synchrotron self-absorption 
(SSA) by formulations given in Li & Chevalier (1999), and 
also free-free absorption (FFA) by formulations in Weiler et al. 
(1986), assuming pre-shocked CSM temperature 7csm = 10 4 K. 
We fix e e = 0.05, and find best-fit parameters of M w and eg to 
the observed radio data (BKC02) by x 2 analysis, as a function 
of the beaming factor b. The radio flux is showing an evidence 
of modulation, presumably due to the interstellar scattering and 
scintillation (1SS), and the minimum y 2 is unacceptably large 
without this effect taken into account. Here we calculate ISS 
modulation index with parameters given in BKC02, and it is 
added to the observational flux errors as quadratic sum. The 
wind velocity is fixed to V w = 10 3 km/s, and different values of 
V w simply rescale M w via the CSM density (oc M w /V w ). The 
change of e e is also mostly canceled by scaling of M w , except 
for the strength of the free-free absorption. We checked that 
changing e e by one order of magnitude does not affect conclu- 
sions derived below. The magnetic field strength can be ex- 
pressed as: B = 0.21 M^V'^WO.Ol) 1 / 2 G. 

The best-fit M w and e B , as well as the y 2 value, are given 
in Fig. 2. Here we used two extreme values of -f m : a low 
value -f m j = 7 and a high value j m j, = 1 + fi e m p (j- l)/m e for 
the thick and thin lines, respectively. The former corresponds 
to a case that the electron minimum energy simply reflects the 
velocity of the shock, while the latter to a case that the kinetic 
energy of ions is efficiently transfered to electrons. We also 
used three values of p = 2.2,2.5, and 2.8. The characteristic 
synchrotron frequency (v m ) corresponding to j m , the SSA fre- 
quency (iAj Sa ), and the FFA frequency (^ff a ) at the day 7 in these 
results are given in Fig. 3, but only for the p = 2.2 case. The 
y 2 degree of freedom is «dof = 24-3 = 21, and the minimum 
reduced y 2 = y 2 /«dof is less than the unity, i.e., an acceptable 
fit. The confidence limit projected on the parameter b can be es- 
timated by a region where Ay 2 , i.e., difference of y 2 from the 
minimum, is smaller than a certain value; Ay 2 < 0.19 and 0.32 
for 95.4 and 99% C.L., respectively, assuming a pure Gaus- 
sian statistics (e.g., Press et al. 1992). Therefore, we conclude 
that a mild beaming b ~ 0.1 is marginally allowed and stronger 
beaming is excluded for the possibility (i). The flux evolution 
and comparison with observed data are shown in Fig. 4, for the 
best-fit models with b = 0.1 and 1. The result in the isotropic 
case (b = 1) is similar to that of BKC02, as it should be. 

A general trend seen in Fig. 2 can be understood as fol- 
lows. When jet is more strongly collimated, the amount of 
CSM swept-up by the jet becomes smaller, and hence higher 
mass loss rate is required to compensate this. However, the ob- 
served spectral feature is mostly explained by SSA, and hence 
magnetic field must become smaller to keep SSA frequency at 
the observed value. This explains behaviors between b = 0.1- 
1 . However, FFA becomes significant when M w becomes very 
large at b < 0.1. The observed data are not fitted well only by 
spectral break by FFA, because the early rise of radio flux due 
to decreasing optical depth is more rapid than SSA (see Weiler 
et al. 1986 for radio supernovae showing this feature), and it 
does not fit the observed slow rise of radio flux at 1.43 GHz. 
As a result, M w cannot increase significantly with decreasing 
b at b < 0.1, and SSA frequency is always higher than FFA 
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for the best fit models (see Fig. 3). Because of this constraint, 
the x 2 rapidly increase with decreasing b smaller than ^0.1. 
The characteristic frequency v m is much smaller than GHz for 
the low 7 m value, but it becomes close to GHz for the high j m 
value. This provides another freedom to the modeling of the 
observed spectrum, and hence slightly better fits. 

4.2. Future emission by interaction with wind bubble or ISM 

When the jet or fast ejecta is interacting with CSM of stellar 
wind profile, the synchrotron flux decreases with time even if 
it is not decelerated. However, it should eventually enter to a 
region where the density is rather uniform. It is either interstel- 
lar medium (ISM) that has not been affected by the progenitor 
star, or uniform density region of CSM such as the stellar wind 
bubble composed by the shocked wind between the contact dis- 
continuity with ISM and wind termination shock (Weaver et 
al. 1977; Koo & McKee 1992). If the jet is not yet decel- 
erated in such uniform density region, the synchrotron radia- 
tion flux should be simply proportional to the mass of swept-up 
CSM/ISM matter, and hence should rapidly increase with time 
as oc t 3 , until the jet sweeps up CSM/ISM matter of a compa- 
rable mass with the jet itself and deceleration starts. 4 (Note 
that the cooling frequency of electrons is much higher than 
the radio bands, and hence cooling does not affect the radio 
flux.) In the jet hypothesis considered here, this occurs at a time 

fdec ~ 13.9/7_| //3 (n e xt/cm" 3 )" 1 ' /3 yr after the explosion. Here we 
assumed that the uniform density region is composed mostly by 
hydrogen, and n ext is the hydrogen number density. Therefore, 
once the jet enters to the uniform external medium, we expect 
to see a rapid increase of the radio flux in a time scale of years. 

The predicted radio flux by this process is shown in Fig. 5, 
where we used b = 0.03 (jet opening angular radius of <~ 10°), 
p = 2.2, e e = 0.05, and e B = 0.01, which are reasonable values for 
GRB afterglows (Panaitescu & Kumar 2001), as well as radio 
supernovae. We suppose the possibility (ii), and we adopted 
M w = 5 x 10" 7 M Q /yr, which is found by BKC02 by isotropic 
modeling. The calculation is stopped at t = fd ec , beyond which 
deceleration must be taken into account. The early radio flux 
by interaction between the jet and wind-profile CSM is smaller 
than observed with these parameters, as required for the possi- 
bility (ii). Here we assumed that the density profile becomes 
uniform beyond a radius r ext = 10 17 or 10 18 cm, and used three 
values of n ext (0.1, 1, and 10 cm -3 as typical values for ISM as 
well as found in GRB afterglows). There is a discontinuity in 
the density at r ext depending on M w and n ext , but such discon- 
tinuity is also expected in reality for stellar wind bubbles (e.g., 
Weaver et al. 1977). 

The transition radius r ext is not easy to estimate, but a cal- 
culation of CSM density profile of presupernova Wolf-Rayet 
stars indicates r ext ~ 10 cm (Ramirez-Ruiz et al. 2001). Ob- 
servations of the prototype ring nebular NGC 6888 harboring 
a WR star indicate that the wind termination shock is at < 
1-3 pc within the hot stellar wind bubble (Wrigge, Wendker, 
& Wisotzki 1994; Moore, Hester, & Scowen 2000). On the 
other hand, the radio light curves of SN 1998bw can be fit- 
ted with the wind density profile up to ~ a few x 10 17 cm (Li 
& Chevalier 1999). Another indirect estimate of r ext is pos- 
sible from GRB afterglow light curves, assuming that the en- 
vironment of SN 2002ap is similar to that of typical GRBs. 

4 It should be noted that the flux does not increase in the case of stellar wind 
electrons increases as oc t. This is because the magnetic field and hence syncl 
energy density and shocked matter. 



At least some of GRB afterglow light curves can be fit better 
by a model with homogeneous ambient density profile rather 
than the wind profile (Chevalier & Li 2000; Panaitescu & Ku- 
mar 2002), and such afterglow light curve seems to start at 
about 0.1-1 day after the burst, as seen in the latest GRBs 
(021004, 021211) with very early optical detections (Holland 
et al. 2003; Li et al. 2003; Fox et al. 2003). The ob- 
server's time is related to the location of the shock of afterglows 
as r= 3.0 x 10 17 (£ iso /10 52 erg) 1 / 4 (n ext /lc m - 3 r 1 / 4 (t obs /0.1 d) 1/4 
cm (e.g., Blandford & McKee 1976), where E\ so is the isotropic 
equivalent energy of GRBs, and beyond this distance the den- 
sity profile around GRBs seems uniform (see also Holland et al. 
2003). Therefore, a plausible scale of r ext is ~ 10 17 -10 18 cm. 

It should also be noted that the speed of isotropic supernova 
ejecta producing the observed radio flux is inferred to be ~ 0.3c 
(BKC02), which is comparable or higher than the inferred jet 
speed of 0.23c. Then the isotropic ejecta producing the ob- 
served radio flux may be propagating faster than the jet. In this 
case the jet is not directly interacting with CSM and early ra- 
dio flux prediction in the unshocked stellar wind region is not 
valid. However, the isotropic supernova ejecta with a speed 
higher than 0.3c is probably a tiny part of the total kinetic en- 
ergy of the supernova, and hence it will start to be decelerated 
much earlier than the jet. Then eventually the jet will become 
the fastest and most distant component among the supernova 
ejecta, and the radio flux prediction in Fig. 5 is valid in such a 
later epoch. 

As expected, the radio flux shows a rapid increase at around 
tat = W(/30 = 4.6(r ext /10 18 cmX/3/0.23)- 1 yr. The radio flux 
would become as strong as ^0.1 Jy. This is not surprising 
when it is compared with radio GRB afterglows, which have 
a similar energy scale but much higher expanding speed; typ- 
ical GRB afterglows would have radio flux of <~ Jy for years 
after the burst, if it is placed at a distance of SN 2002ap (Totani 
& Panaitescu 2001). The expected radio flux for SN 2002ap 
in the near future is strong enough to detect by a long-term 
monitoring. Furthermore, jets expanding to two opposite direc- 
tions (both having 6* t, s ~ 90°) will have an angular separation 
of # S ep = 20(f/5 yr)(/3/0.23) mas, which is easy to spatially re- 
solve by VLBI observations, as proven for SN 1993 J (Bartel et 
al. 1994;Marcaideetal. 1995,1997). The jet directions should 
be perpendicular to the observed polarization angle, giving a 
decisive proof for the jet hypothesis if it is actually detected. 
On the other hand, non-detection in next tens of years would 
not necessarily exclude the jet hypothesis; if the density of the 
uniform region is low (n ext < 0. 1 cm -3 ) and r ext is large, it would 
take a long time (> a few tens of years) until the radio emission 
becomes detectable again. 

5. DISCUSSION 

5.1. Implications for the GRB-SN Connection 

The jet inferred from SN 2002ap, having a similar energy 
scale to GRB jets, should still be regarded as a hypothesis, and 
we should not over interpret it at this time. However, since we 
have shown that the jet is physically possible and it can be tested 
in the future, it is interesting to think what the jet would mean 
and how it would fit to the other observational facts obtained 
so far concerning the GRB-SN connection. Here we suggest a 
plausible picture for the GRB-SN connection putting together 

:rnal profile, as seen in the previous section, even if the total energy of shocked 
•on emissivity decrease with time by the scaling assumed between the magnetic 
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the jet from SN 2002ap and other observational results, assum- 
ing that the jet of SN 2002ap is real. 

The similarity of the jet energy scale indicates that the SN 
2002ap jet is produced by the same mechanism as that for 
the cosmologically distant GRBs. Then it may be called as 
a "failed" GRB, or another possibility is an off-axis GRB, as 
mentioned in § 1 . Fortunately, we can reject the latter possi- 
bility by radio observations made so far. If an off-axis GRB 
is the case, then we expect even earlier re-appearance of ra- 
dio emission than discussed in §4.2, which should show even 
faster expansion with a velocity of ~ c, than the jet having 
a velocity of ~ 0.23c. In fact, it is exactly what is called as 
"orphan afterglows": a GRB whose jet direction is away from 
an observer, which can be seen only by late time, less colli- 
mated afterglow emission. [For discussions on detectability of 
radio emission from nearby GRB remnants, see also Paczyhski 
(2001).] As shown in Totani & Panaitescu (2002), typical GRB 
afterglows viewed at a large angle from the jet center at day 
~100 should have flux of - 10-100 Jy, V ~ 16-18, and - (1- 
10) x 10~ 12 erg cm" 2 s _1 , in radio (5GHz), optical, and X-ray (1 
keV) bands, respectively, at the distance to SN 2002ap. These 
fluxes are too faint to detect if it is placed at cosmological dis- 
tances, but thanks to the close distance to SN 2002ap, we should 
easily detect an orphan afterglow of SN 2002ap. Then, the lat- 
est radio data of SN 2002ap about 200 days after the explosion, 
that is still at ~ 0.1 mJy level showing no evidence of flaring 
up (Berger et al. 2003), already exclude the off-axis GRB pos- 
sibility. 

As discussed in §3.4, decay gamma-ray of 56 Ni is the most 
likely ionization process of the jet from SN 2002ap, but the 
required amount of 56 Ni in the jet and/or outside the photo- 
sphere is difficult to explain only by hydrodynamical instability. 
Hence, this indicates that the jet is formed and ejected from the 
central region of the core collapse. On the other hand, we do not 
expect sufficient amount of 56 Ni for ionization, if the jets are 
formed at outer envelope of the star, such as trans-relativistic 
acceleration of shock wave when it passes through the steep 
density gradient of stellar surface, as studied by Matzner & Mc- 
Kee (1999) and Tan, Matzner, & McKee (2001). 5 

On the other hand, considering the similarity between SN 
2002ap and 1998bw / GRB 980425, we may also want to estab- 
lish a consistent picture including GRB 980425 which is pecu- 
liar in a few points compared with GRBs found at cosmological 
distances. GRB 980425 has soft spectrum and smooth tempo- 
ral profile (such as long spectral lag), and its total isotropic- 
equivalent energy is about 10 4 times smaller than typical GRBs 
at cosmological distances (Bloom et al. 1998; Norris et al. 
2000). One possible interpretation of these properties is that 
this is a typical GRB but we observed this event at a slightly 
off-axis direction from the jet center (Nakamura 1999; Ioka & 
Nakamura 2001; Salmonson 2001; Granot et al. 2002). How- 
ever, if GRB 980425 is a normal but off-axis luminous GRB, 
we expect orphan afterglow radiation after it became mildly- 
relativistic. Although optical afterglow radiation may be hid- 
den by the brightness of SN 1998bw (Granot et al. 2002), radio 
afterglow flux expected for typical GRBs at a distance to SN 
1998bw is 0.4^1 Jy at t ~ 100 days (Totani & Panaitescu 2002), 
which is still much larger than actually observed (~ 10 mJy 
at day 100), even though SN 1998bw was much brighter than 
SN 2002ap. In addition to this, Norris (2002) identified a sub- 



class of GRBs including GRB 980425, whose members have 
low luminosity, long spectral lags, and soft spectrum. These 
low-luminosity GRBs tend to have just a few wide pulses while 
nearly all high-luminosity GRBs have many, narrow pulses, in- 
dicating that only the viewing angle effect cannot explain the 
difference between the low- and high-luminosity GRBs. 

Here we note our result that the radio emission from SN 
2002ap cannot be strongly collimated. Polarization was not de- 
tected for the strong radio emission of SN 1998bw, indicating 
that the radio emission is produced by isotropic ejecta (Kulkarni 
et al. 1998). The high velocity shell responsible for radio flux 
of SN 1998bw shows an evidence of deceleration (Li & Cheva- 
lier 1999), and hence the kinetic energy of the shell should not 
be much larger than ~ 10 49 erg. GRB 980425 associated with 
SN 1998bw has isotropic equivalent energy of about ~ 10 48 erg, 
which is very different from the standard energy scale of cos- 
mologically distant GRBs, but is close to that of isotropic radio 
emitting shell. 

Tan et al. (2001) have shown that the peculiar GRB 980425 
can be produced by the isotropic fast shell produced by the hy- 
drodynamical shock acceleration, which is the same shell pro- 
ducing radio emission. The smooth temporal profile of GRB 
980425 may also be expected by isotropic shock acceleration 
process. Even if explosion is highly asymmetric at the center, 
a two-dimensional simulation by Maeda et al. (2002) indicates 
that the outer low-density region is rather isotropic. The ef- 
ficiency of the shock acceleration sensitively depends on the 
outer density profile of stellar envelope (Tan et al. 2001), and 
it may not be surprising even if SN 1998bw and SN 2002ap 
produced very different velocities and energies of fast moving 
isotropic ejecta. 

Then, the apparent discreteness between GRB 980425 and 
other GRBs may be understood by two distinct driving pro- 
cesses of GRBs: cosmologically distant GRBs are strongly col- 
limated jets with an energy scale of ~ 10 50-51 erg, which is 
produced along with ejection of 56 Ni outside the photosphere 
by the central activity of core collapses (refereed as type 1 for 
convenience), while GRB 980425 is produced by much less en- 
ergetic, isotropic fast ejecta, which could be produced by hy- 
drodynamic shock acceleration at outer layer of exploding stars 
(type 2). There is a possibility that SN 1998bw was successful 
also as off-axis type 1 GRB, but this possibility is not favored 
by the no detection of radio orphan afterglow showing an en- 
ergy scale of ~ 10 51 erg for ~ 300 days, as mentioned above. 

Other past events of type Ic supernovae may also have pro- 
duced energetic jets with a velocity of ~ 0.2c, like SN 2002ap. 
Such jets, if exist, may now be emitting strong radio emission 
after sweeping up enough amount of ISM/CSM, showing jet- 
like morphology. We note that the time scale of emergence in 
radio wavelength could be much longer for the case of a failed 
GRB with low-speed jets, than considered in previous publi- 
cations for off-axis GRBs or GRB remnants (Paczyhski 2001). 
On the other hand, it would be shorter than the time scale of es- 
tablishment of normal radio supernova remnants, correspond- 
ing to the transition from the free expansion to Sedov-Taylor 
phase (~ 100 yr). Re-examination and new follow-up of such 
past type Ic events are encouraged. 

The picture of GRB-SN connection presented here predicts 
that all type 1 GRBs should be associated with energetic type 
Voile supernovae (as confirmed by SN 2003dh / GRB 980329 



5 Following Tan et al. (2001), we use "hydrodynamical shock acceleration" for this phenomenon, which must not be confused with the Fermi acceleration of 
cosmic-ray particles. 
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after the submission of this paper), but we expect a variety of 
supernova luminosity that may not be correlated with GRB lu- 
minosity or the jet energy. The supernova luminosity is de- 
termined by the amount of synthesized 56 Ni within the pho- 
tosphere where the radioactive energy is converted into optical 
photons. As we suggested, the jet formation activity may eject a 
significant amount of 56 Ni well outside the photosphere, where 
56 Ni cannot contribute to the optical luminosity. This is not in- 
consistent with the results of search for supernova evidence in 
GRB afterglows (Bloom et al. 1999; Galama et al. 2000; Re- 
ichart 2001; Bloom et al. 2002; Garnavich et al. 2003; Price et 
al. 2003). 

5.2. On the X-ray line features in GRB afterglows 

Emission line features of iron (or nickel) and multiple-alpha 
nuclei (Mg, Si, S, Ar, Ca, etc.) are often found in X-ray GRB 
afterglows on a time scale of ~ a day (Piro et al. 1999, 2000; 
Yoshida et al. 1999; Antonelli et al. 2000; Reeves et al. 2002; 
Butler et al. 2003). Theoretical explanations mostly fall into the 
two categories: (1) geometry-dominated (GD) scenario where 
the time scale of ~ a day is attributed to the photon propaga- 
tion time, and this scenario needs a supernova explosion that 
occurred weeks prior to a GRB (the supranova model, Vietri 
et al. 2001), dense preburst circumstellar material (Weth et 
al. 2000; Kotake & Nagataki 2001) or a distant reflector of 
GRB/afterglow emission such as a e ± -pair screen (Kumar & 
Narayan 2003); and (2) engine-dominated (ED) model where a 
long-lived energy source left over in the center of the star af- 
ter the GRB is ionizing the matter around it (Rees & Meszaros 
2001). Most of these explanations assume that ionization pro- 
cess is photoionization by GRB/afterglow radiation or long- 
lived remnants, while an alternative is shock heating around the 
center (Bottcher 2000). However, all these scenarios have one 
or more problems (Lazzati, Ramrez-Ruiz, & Rees 2002; Ku- 
mar & Narayan 2003), and other possible explanations are still 
worth seeking for. 

Here we propose that the ionization by radioactivity of de- 
caying 56 Ni in mildly optically thin region, which has been sug- 
gested as the ionization mechanism for the SN 2002ap jet, can 
be considered as a new explanation for the X-ray line features 
of afterglows. First let us check the energetics. The luminos- 
ity of decaying energy is given by ~ 1 .0 x 10 44 (M N i/M Q ) erg/s 
within ~ = 8.1d after the explosion. We need several solar 
mass of 56 Ni to explain the luminosities of X-ray line features; 
this is admittedly large compared with typical supernovae, and 
here we have assumed 100% efficiency for the energy conver- 
sion from radioactive decay into line photons, which seems 
rather extreme. However, we know that SN 1998bw produced 
about ~ IMq of nickel, and it seems not very unlikely that more 
massive and energetic supernova (or hypernova) events produce 
even more nickel than SN 1998bw and eject it out to relatively 
optically thin region, via the process of jet formation for type 1 
GRBs. If radioactive gamma-rays lose their energy dominantly 
by ionizing neaby ions, efficiency of line production could be 
close to the unity in some preferred situations. If all 56 Ni ions 
are fully ionized, the radioactive decay by electron capture is 
impossible (McLaughlin & Wijers 2002), but repeated decays 
and ionizations would be possible if recombination rate is suf- 
ficiently high. This is a complicated process, and clearly more 
careful and quantitative study is required to support this specu- 
lation. 

This scenario can be considered as a new type of ED expla- 



nations that do not need a supernova prior to a GRB, but has one 
important difference from previous ones; this scenario does not 
require strong radiation source in X-ray band for ionization, but 
high energy electrons produced by scattering of decay gamma- 
rays directly ionize ions, and hence we might expect high equiv- 
alent width with weak continuum radiation. In this way the ma- 
jor problem of ED scenarios, i.e., that we should have directly 
observed stronger ionizing radiation or bremsstrahlung contin- 
uum of shock-heated matter than actually observed (Lazzati et 
al. 2002; Kumar & Narayan 2003), might be avoided. 

The time scale of line production should be determined by 
evolution of physical conditions by expansion, such as optical 
depth, density, and recombination rates. Note again that ra- 
dioactive ionization should be efficient only in a small range of 
optical depth where it is mildly optically thin, and then we ex- 
pect a shorter time scale of line production than typical peak of 
supernova light curve that is determined by diffusion of radia- 
tion within the photosphere, and it might become even shorter 
than the time scale of 56 Ni decay. 

6. SUMMARY AND CONCLUSIONS 

We have shown that the jet hypothesis proposed by Kawabata 
et al. (2002) for SN 2002ap based on their spectropolarimetric 
observation at day ~ 13 is physically possible and consistent 
with all observations. We estimated the jet mass (<~ 0.01M Q ), 
direction (approximately perpendicular to the observer), and 
energy (~ 5 x 10 50 erg) by a fully relativistic treatment. The 
large jet energy does not contradict with the weak radio flux, 
since the jet is almost freely expanding and jet material is not 
yet shocked. The total energy of electrons inferred from syn- 
chrotron radio flux, ~ 1 45 erg, only reflects that of CSM swept- 
up and shocked by jet or high-speed ejecta from SN 2002ap. 
The jet becomes optically thin on a time scale of 5-10 days, 
and weak continuum polarization found by the earlier VLT ob- 
servation at day ~ 5 can be explained by high optical thickness 
of the jet to electron scattering. 

The jet must be substantially ionized to reproduce the red- 
shifted polarized continuum. Rapid adiabatic loss of the jet 
internal energy should have cooled down the jet material be- 
low the temperature required for the jet to be kept ionized, and 
hence external heating or photoionization source is necessary. 
We have shown that photoionization is quite unlikely by any 
sources, and the most likely ionization process is heating by 
decaying 56 Ni when the jet is mildly optically thin. Still, ion- 
ization is not easy; the jet must have a mass fraction of order 
unity of 56 Ni, or a larger amount of 56 Ni outside the photo- 
sphere that is at least comparable with that within the photo- 
sphere (~ 0.07M Q ), which is difficult to explain simply by hy- 
drostatic instability. This result indicates that the jet must be 
formed and ejected from the central region of the core collapse. 

We examined whether the observed radio emission can be 
explained by CSM swept-up by the jets, and we found that the 
radio data favor isotropic ejecta. Jets with b = (£lj e t/47r) ^0.1 
is marginally possible to explain the radio data, and stronger 
collimation is excluded as the origin of the observed radio emis- 
sion. If the jet is strongly collimated, the radio emission must 
be from isotropic fast ejecta that is a different component from 
the jet. 

The jet should be regarded still a hypothesis, but we pre- 
dict that, if the jet hypothesis is the case, the jet would even- 
tually pass through the unshocked presupernova wind region 
(p oc r" 2 ), and enter to a region where the CSM/ISM density is 
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rather uniform. The jet mass is large enough not to be deceler- 
ated until it sweeps up a comparable mass of ISM/CSM in a few 
to ~10 years. Then we expect a rapid increase of radio flux that 
should easily be detected by future long-term radio monitoring, 
with reasonable parameters for ISM/CSM density and profile. 
Fortunately, such radio emission can be resolved spatially by a 
VLBI observation, and the morphology and jet direction rela- 
tive to the observed polarization angle would give a clear proof 
of the jet hypothesis. 

Although it is a speculation that would become effective 
when the jet hypothesis is confirmed by future observations, we 
discussed how the inferred jet would fit to the our knowledge 
about the GRB-SN connection, taking into account various ob- 
servational facts obtained in the past. To explain all the obser- 
vations consistently, we suggest two distinct classes of GRBs 
by different formation processes but from similar core collapse 
events; cosmologically distant GRBs are produced by strongly 
collimated jets having an energy scale of ~ 10 50 " 51 erg, which 



are produced and ejected from the central region of the core 
collapse. SN 2002ap jet can be considered as a failed GRB of 
this class, with much larger baryon load than typical successful 
GRBs. 56 Ni is ejected from the center along with the jet, which 
is responsible for the ionization of the SN 2002ap jet. Such 56 Ni 
may also explain the X-ray line features often found in GRB 
afterglows, depending on the amount of 56 Ni and physical con- 
ditions, although more quantitative study is required to verify 
this possibility. On the other hand, the peculiar GRB 980425 
and radio emission from SN 1998bw and SN 2002ap that seem 
isotropic may be produced by isotropic and much less energetic 
ejecta, which could be formed by the hydrodynamical shock 
acceleration at the surface of an exploding star. 
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FIG. 1 . — Top and middle panels: the jet velocity f}, Lorentz factor 7, the jet mass Mj ct , and the jet kinetic energy £j et required to explain the observed redshift 
(z = 0.3) and polarization degree (fp = 1.8 X 10~ 3 ), as a function of the observer's direction obs measured from the jet direction. Bottom panel: the polarization 
degree (/p, p) and redshift (z op ) expected for the opposite jet having the direction of (180° - obs ), with the same jet parameters shown in the top and middle panels. 
Observed redshift and fp are marked by horizontal lines. 
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b = (Q. t /47T) 



FIG. 2. — Presupernova mass loss rate M w (top panel), fractional energy density of magnetic field (middle panel), and reduced \ 2 obtained by \ 2 fitting of the 
radio emission expected from swept-up CSM by the jet to the observed data. The solid, dashed, and dot-dashed curves are for different electron power index as p = 
2.2, 2.5, and 2.8, respectively. The low and high values of the minimum electron Lorentz factor, (7 m j and 7„j,), are used for thick and thin curves, respectively. 
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FIG. 3. — The characteristic synchrotron frequency v m (dashed), synchrotron self absorption frequency i/ ssa (solid), and free-free absorption frequency ug a (dot- 
dashed) at the day 7, for the best-fit models to the observed radio flux and spectrum of SN 2002ap. The electron power index is fixed to p = 2.2, and thick and thin 
curves are for the low and high values of 7„: y m j and "/ m h , respectively. 
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FIG. 4. — The radio flux evolution of the best-fit models for (b./y m ) = (l,7 m .;), (1,7m.*), (0.1, 7 m .;), and (0.1, 7 m A ), for the solid, dotted, dashed, and dot-dashed 
lines, respectively. The observed data are from BKC02, with the thick error bars showing observational errors while the thin error bars include the ISS modulation 
by quadratic sum to the observational errors. 
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FIG. 5. — Prediction of radio flux by CSM/1SM swept up by the jet in the future. Here we assumed that the observed radio flux (data points from BKC02, 
error bars not including ISS modulation) is from CSM swept up by isotropic supernova ejecta that is a different component from the jet. The parameters of 
(n e xt/cirr 3 ,rexi/cm) = (1, 10 17 ), (0.1, 10 17 ), (10, 10 17 ), and (1, 10 ) are used for the solid, dotted, dashed, and dot-dashed lines, respectively. 



